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A SYMBOLIC TREATMENT OF RIORDAN ARRAYS 


jose agapito, Angela mestre, pasquale petrullo, and maria m. Torres 


Abstract. We approach Riordan arrays and their generalizations via umbral symbolic methods. This 
new approach allows us to derive fundamental aspects of the theory of Riordan arrays as immediate 
consequences of the umbral version of the classical Abel’s identity for polynomials. In particular, we 
obtain a novel non-recursive formula for Riordan arrays and derive, from this new formula, some known 
recurrences and a new recurrence relation for Riordan arrays. 

1. Introduction 

A Riordan array is an infinite lower triangular matrix R = {r nt k)n,k£K ordinarily described by a pair of 
generating functions (d(t), h(t)) such that rf(O) 7 ^ 0, h(0 ) = 0, h'(0 ) 7 ^ 0, and r n ,k = [t n } d(t) h{t) k , where 
[t n ] is the operator which gives the n th coefficient in the series development of a generating function. If 
d(t) = X^ra>o T> h(t) = E n >i h n ^j, and r n ^ = n\ [ t n ] d(t) , then R is called exponential Riordan 
array. More generally, given a sequence (c n ) of nonzero numbers and taking d(t) = XE>o hit) = 

XE>i h n t—, and r n ^k = c n [t n ]d{t) h ^ , R is called a generalized Riordan array with respect to (c n ) 
[551 . Thus, for example, the Riordan array whose entries are the binomial numbers (£) can be described 
by the pair (y^y, jzii) (c„ = 1 : ordinary presentation) or by (e 4 ,f) {c n = n!: exponential presentation). 
Riordan arrays form a group under matrix multiplication. The literature about Riordan arrays is vast 
and still growing and the applications cover a wide range of subjects, such as enumerative combinatorics, 
combinatorial sums, recurrence relations and computer science, among other topics u ns na na gs 
HU HU [SH [321 [Ml M- Formally, ordinary Riordan arrays are a formulation of the 1-umbral calculus, 
whereas exponential Riordan arrays are a formulation of the n!-umbral calculus. In fact, the classical 
umbral calculus, as given in m consists of a systematic study of a certain class of univariate polynomial 
families ( Sheffer sequences ) by employing linear operators on polynomials. An extension of the classical 
umbral calculus to infinitely many variables can be found in [5]. Sheffer sequences, provided with a 
particular non-commutative multiplication ( umbral composition), form a group which is isomorphic to 
the Riordan group [12,. 

The purpose of this paper is to give a promising new symbolic treatment of Riordan arrays, based 
on a renewed approach to umbral calculus initiated by Rota and Taylor in [30,. This renewed approach 
makes no use of operator theory. The symbolic techniques developed from this new approach have been 

fruitfully applied to a wide range of topics [7[ ISJ121 [TUI HU 1201122112211211122112U12U [Ml EH] • In particular, 

Di Nardo, Niederhausen and Senato [7] gave a representation of exponential Riordan arrays as a result of 
their symbolic handling of Sheffer polynomials. Definition 13.11 given later on this paper is a normalized 
version of their formulation. A treatment of Riordan arrays under the more general context of recursive 
matrices using the more established approach to umbral calculus m can be found in mm- 
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A key tool in our symbolic approach to Riordan arrays is the urnbral version (formula ©) of the 
following binomial Abel identity for polynomials 

(i) (* + y) n = p Q {y + ka) n ~ k x(x - kaf- 1 . 

The urnbral analogue of © was first stated in [9] (without a proof) and later obtained from an urnbral 
version of the Lagrange inversion formula in [7]. An elementary proof of the urnbral version of Abel’s 
identity has been recently given by Petrullo [23; :24j . 

The main contribution of this paper is to show that well-known properties of all Riordan arrays are 
easily derived from the urnbral version of Abel’s identity. This is a nontrivial aspect from a theoretical 
point of view. It is not immediate by means of the traditional methods for dealing with Riordan arrays. 
It is worth mentioning that, in the classical context of generating functions, Sprugnoli [34] had already 
used Riordan arrays to obtain several combinatorial formulas that generalize the classical Abel’s identity. 
Therefore, our point of view is conceptually different from Sprugnoli’s approach. 

The paper is organized as follows. Section [2] recalls the basics of Rota and Taylor’s classical urnbral 
calculus m- Along the way, two equivalent urnbral versions of Abel’s identity for polynomials and an 
urnbral version of the classical Lagrange inversion formula are given. Section [3] states the fundamental 
aspects of the theory of Riordan arrays (of exponential type) in urnbral terms. Also, connections with 
Sheffer sequences are recast in the new urnbral symbolic setting. Furthermore, some recurrence relations 
of Riordan arrays are stated. Section [4] tests our urnbral approach to Riordan arrays on some concrete 
classical examples, providing both urnbral and traditional formulas. Section [5] extends the discussion 
given in Section [3] to w-Riordan arrays, where w is any umbra with nonzero moments. In this generalized 
context, we obtain an important urnbral formula (Theorem l5.3l) as a direct consequence of the urnbral Abel 
identity ©. One of our main results is a novel non-recursive formula (Theorem 15.41) . which is a direct 
application of Theorem 15.31 Known recurrence relations and a new recurrence formula for w-Riordan 
arrays (Theorems 15.5115.61 and 15.711 are derived from Theorem 15.41 

2. Preliminaries 

This section starts with a brief review of the essentials of the classical urnbral calculus that underlies 
much of this paper. We recall the intimate relation that exists between the urnbral version of Abel’s 
identity © and the Lagrange inversion formula. This connection is not clear at first sight from the 
traditional point of view using generating functions (see for instance [ 6 ] for the common view of these 
two classical formulas). The urnbral versions of the binomial Abel identity and of the Lagrange inversion 
formula play an important role in our symbolic approach to Riordan arrays (see Sections [3] and 0 . 

2.1. The basics of classical urnbral calculus. To perform classical urnbral calculus we need to have: 

(i) A commutative integral domain R with identity 1 (assume R = C[t], x possibly multivariate). 

(ii) An alphabet A = {a,/ 3 , 7,...} with enough letters called umbrae. 

(iii) An iJ-linear functional E: i?[A] — > R such that E[ 1] = 1 and E[a 1 ft ■ ■ ■ y fc ] = E[a l ]E[^] ■ ■ ■ E[y k ], 
for pairwise distinct umbrae a,/3,... ,7 (uncorrelation property). The map E is called evaluation. 

I? [A] simply denotes the polynomial ring in the set A of commuting indeterminates, with coefficients in R, 
endowed with formal addition and multiplication. An element p £ i?.[A] is called urnbral polynomial. The 
support of an urnbral polynomial p £ I? [A] is defined as the set of all umbrae that show up as variables 
of p. Two urnbral polynomials p, q £ R[A\ are said to be uncorrelated if their supports are disjoint. 
In particular, distinct umbrae are always uncorrelated. Ordinary equality defines a trivial equivalence 
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relation on R\A\. Additionally, two nontrivial equivalence relations can be defined on R[A\. We write 
p ~ q and say that p and q are umbrally equivalent , when E[p\ = E[q\. Moreover, if p n ~ q n for all 
integers n > 0, we say that p and q are similar and write p = q. An umbral polynomial p is said to 
represent a sequence (a n ) in R (with ao = 1) if E[p n ] = a n for all n > 0. In this case we say that a n is the 
n th moment of p and write p n ~ a n . The equivalent relation p = q means that p and q represent the same 
sequence. Often, it is convenient to assume that each monic sequence (a n ) in R is represented by infinitely 
many distinct similar umbrae. In this case we deal with a saturated umbral calculus or, equivalently, with 
a saturated alphabet [301 . We can then rephrase property (ii) in a more precise manner by saying that 
A is saturated. 

Remark 2.1. Although two equal polynomials are similar and two similar polynomials are umbrally 
equivalent, neither converse holds. For instance, consider umbrae a, a', 7 , 7 ' such that a 7 ^ a', a = a' 
and 7 7^ 7', 7 = 7'. Clearly a + 7 / a' + 7', but a + 7 ~ a' + 7'. Moreover, we have a + 7 = a' + 7'. 
Similarly, note that a + a' ~ 2a, whereas a + a' ^ 2a. These examples show that identities in i?[A] 
involving umbral equivalence and similarity encode more information than ordinary equality. 

Given p £ R.[A], we denote by e pt the formal exponential series 

00 j-n 

+ £ R[A][[t]]. 

71=1 

The series above is called the generating function (g.f. for short) of p. We extend the action of E from 
i?[A] to i?[A][[f]] and write f p (t) = E[e pt ] := 1 + E\p n ] ^ £ i?[[i]]. The series f p (t) is called the 

moment generating function (rn.g.f. for short) of p. As done for umbral polynomials, we may also define 
umbral equivalence on 7?[A][[f]] and write e pt ~ e qt if and only if E[e pt ] ~ E[e qt ], so that e pt ~ 

It is not difficult to see that f p (t) = f q (t) if and only if p = q. We will often use the terms rn.g.f. 
and g.f without distinction throughout the paper. For later reference, we summarize in Table [T| useful 
information on some fundamental umbrae. 


Umbra 

a. 

f a (t) = E[e at ] 

E[a n ] = n\[t n }f a {t) 

Augmentation 

£ 

1 

1,0,0,... 

Unity 

V 

e t 

1,1,1,... 

Singleton 

X 

1 + t 

1,1,0,... 

Bell 

(3 

e 6 * -1 

1,1, B 2 , B 3 ,... (Bell numbers) 

Bernoulli 

L 

t 

e t — 1 

1, b±, b 2 , ■ ■ ■ (Bernoulli numbers) 

Boolean unity 

V 

1 

1 -t 

0!,1!,2!,3!,... 


Table 1. Some fundamental umbrae. 


On the one hand, note that f a+1 {t) = I?[e^“ +7 ^] = .Efe^e 74 ] = E[e at ] A[e 7 *] = / a (t)/ 7 (t) if a 7 ^ 7 
(uncorrelation property). On the other hand, clearly / 2 a (f) = E[e^ 2a ' ,t ] = E[e a ^ 2t ' ) ] = f a (2t) 7 ^ (/ a (t)) 2 . 
However, if a 7 ^ a' and a = a', then f a + a ' (t) = ( f a (t )) 2 . The inequality f a + a '(t) 7 ^ ha(t) reflects the 
fact that the umbral polynomials a + a' and 2a are not similar, as was pointed out in Remark 12.11 We 
stress this difference and set 2.a = a + a' (with a/a' and a = a'). Thus, for all n > 0, we have 

( 2 .ar * (a + aT * £ (a')"" fc ^ £ Qa fe a„_ fe . 
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Up to similarity, the auxiliary umbra 2 .a is well defined. In general, given x £ RU A, we shall implicitly 
define the auxiliary umbra x.a by means of its g.f.; eU-“)* ~ that is, ( x.a) n ~ n\[t n ](f a (t)) x . 

In particular, for any 7 , a £ A, we have e^ 7 ‘“^ ~ (/ Q (t )) 7 — e 7 ( los -GW); namely, / 7 . a (f) = / 7 (log/ a (t)). 

Remark 2.2. Let a™ ~ a n for all n £ IN. The moments of 7 .a are explicitly given by the formula 

n 

( 7 .a) n ~ • • •, a n -i+ 1 ), 

i =0 

where ( 7 ); = 7(7 — 1) • • • (7 — i + 1) and the < 22 ,..., a n —i+ 1 ) are the partial Bell exponential 

polynomials (see [9] for details). Note that 7 .a ~ 7 a but 7 .a ^ 70 . 

The dot product of umbrae is in general non-commutative. Nevertheless, it satisfies several useful 
identities which can be easily verified via generating functions. 

Example 2.1. Given a € A, r £ R and the umbrae e, j3, v and \ whose g.f.’s are shown in Table [Q we 
have 

a.e = e = e.a , y./J = /3.x = v , ra = a.(rv) , a + e = a = e + a , a + (— l).a = e = (—l).a + a. 

For instance, / s . a (t) = / £ (log/ a (t)) = 1 and / a . e (t) = / Q (log/ e (t)) = / a (0) = 1; hence a.e = e = e.a. 
Similarly, f x .p(t) = f x (\ogf 0 (t)) = f x { loge 6 *- 1 ) = 1 + e* - 1 = e‘ and // 3 . x (t) = fp(}ogf x {t)) = 
f( 3 (log(l + t)) = e e log(1+f) -i = e 1+t_1 = e*; hence y./3 = /3.y = t>. 

By iterating the dot product we can define further auxiliary umbrae. In this way, the umbra 7 ./3.a 
(where /3 is the Bell umbra) is called the composition umbra of 7 with a because it encodes the composition 
of formal series; namely, f lm p. a (t) = fj(fa(t) — l). We can also define new umbrae out of a single umbra. 
For example, a <-1> denotes the auxiliary umbra that satisfies a <_ 1 > ./3.a = \ — a./3.a <_1> . We call a <-1> 
the compositional inverse of a. One can check that / a <-i> (i)— 1 = ( f a (t ) — 1) <-1> . Likewise, we define the 
a-derivative umbra q b as the auxiliary umbra whose moments are given by a™ — na n_1 for all n > 1 . 
It is easy to see that the m.g.f. of a v satisfies f a ^ (t) = 1 + t f a {t). Since (t fait))™ ~ t m for any 

integer m > 0, we have e ( - a ' /3 ' a i>') t ~ e CT I*/<»(*)] ~ J2m =0 Y^kLo am i m - a ) k ^rfr- Comparing the coefficients 
of Uf, we obtain the following binomial-like expansion for the moments of a composition umbra, 

n , 

( 2 ) (<,.p. a „r ^Y.\ k 

k—0 ' 

Remark 2.3. Let f a (t) be the m.g.f. of a. Formula @ is equivalent to 

n 

(3) \t n ]fa(tUt)) = £>*]/„(*) [*""*](/„ (t)) k . 

k =0 

2.2. Umbral Abel identity and the Lagrange inversion formula. Let a, 7 and 0 be any umbrae. 
For all n > 1, we have 

(4) (7 + a) n ~ ^ (7 + k.a) n ~ k a(cr + {—k).a) k 1 (Umbral Abel Identity, Version I). 

k =0 ' ' 

Formula ([4]) is a straightforward (though nontrivial) generalization of the classical Abel identity for 
polynomials in three indeterminates shown in 0 - An elementary proof and further applications of 0 
have been recently given by Petrullo [24]. The umbral Abel polynomials o(a + (—n).a)" 1 are direct 
generalizations of the classical Abel polynomials x(x — na) n ^ 1 for a € R. Such polynomials were first 
studied by Rota, Shen and Taylor [29] in connection with polynomials of binomial type. For the sake of 
clarity, our notation for umbral Abel polynomials slightly differs from the usual convention: er(cr—n.a)” 1 . 


^a k (k.a) n ~ k 
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Umbral Abel polynomials have been applied to cumulants and convolutions in probability theory by Di 
Nardo, Petrullo and Senato jSJ, and further generalized by Petrullo in [221 [23] - For convenience, let us 
now denote by an auxiliary umbra whose moments are given by 

(5) - fJ ( cr + (-n).a)" -1 , n> 1. 

A simple proof of the Lagrange inversion formula immediately follows. 

Theorem 2.1 (Lagrange inversion formula). For any umbrae a, 7 and all integers n > 1, we have 

(6) <r(cr+(-n).a) 71 1 ~ (o./3.a ( ~ 1) ) n . 

Proof. By applying ©, then (|5|), and finally we get 

{K. a -P.Ov)” * J2 (l) a ( a + (- k).a) k ~\k.a) n ~ k ~ (e + o) n ~ o n . 
fc=0 ' ' 

Therefore, = a, which implies that = (J.j3.a { ~ 1) . □ 

A plain translation of © in terms of m.g.f.’s of umbrae gives the traditional Lagrange inversion formula 

( 7 ) [t^jm (JL)" = . 

Now, let rj = M a ,a- Since ig.fi.= o.j3.a l ff 1) .fd.a^ = o , we can rewrite ® as 

( 8 ) (7 + rj.(3.a J> ) n ~ ^ (7 + k.a) n ~ k rj k , n > 0. (Umbral Abel Identity, Version II) 

We shall see in Section [3] that formula © gives an umbral version of the fundamental theorem of Riordan 
arrays. In addition, it follows from © and © the next useful umbral identity, 

( 9 ) (a + = 77^.07,. 

Let us introduce our last umbral tool. By £ 7ia , we denote an auxiliary umbra that satisfies £ 7 , a = — l.A 7iQ! . 
We write := and £„ := £„,„. The umbrae M. a and £ Q play a central role in the umbral theory 
of cumulants © :23j. Since £ Q = —l.a. fi.a { ~ 1) and £ v = y, by © we have y = {a + £ a ./3.a I ,) c = 
(£ a ) T ,.j3.a T) , so that 

( 10 ) (£*L = a^ _1> - 

In particular, we have £fi Q = a. It is for this reason that £ a is named the Lagrange involution of a. 
Now, for all k > 1 denote by S (k) an umbra whose g.f. = 1 + ^j-. By (fTUl) . we have £ (fe) ./3. (£ a )i> — 

&<■*>. Hence, by © and ©, we get (^)(fc.£ a )" _fc ~ (— n.a) n ~ k , so that 

(11) n(k.£ a ) n ~ k - k{-n.a) n ~ k . 

Since ££ a = a and —n.£ Q = n.& a , formula ( 1111 ) is equivalent to 

(12) n(k.a) n ~ k ~ k(n.R a ) n ~ k . 

Remark 2.4. As pointed out by G. Andrews [I], at first glance, the classical umbral calculus seems to be 
merely a convenient notation for dealing with exponential generating functions. However, it reveals its 
power in tackling difficult problems concerning generating functions not easily handled by usual algebraic 
techniques mm- The aforementioned umbral identities (involving dot products, umbral equivalence 
and similarity) show not only more appealing formulas but also, from our perspective, more effective ones 
since they allow for further algebraic manipulations that help us to write in polynomial terms identities 
otherwise stated in terms of the operator [t ra ] and generating functions (compare e.g., © with © and 
© with ©). 
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3 . Exponential Riordan arrays 

In this section we use the umbral syntax developed in Section [2] to encode Riordan arrays. We show 
that all the basic facts about Riordan arrays are easily stated and derived by means of elucidating umbral 
identities. For simplicity, we test the effectiveness of the umbral notation on normalized Riordan arrays; 
namely, those having diagonal entries all equal to 1. We provide a simple umbral characterization of the 
group axioms, describe some important Riordan subgroups, and recover the isomorphism between the 
Riordan and the Sheffer groups. Likewise, recursive properties of Riordan arrays are obtained as direct 
consequences of the umbral Abel identity. These properties will be further generalized in Section [5] 


3.1. The exponential Riordan group. Riordan arrays are a generalization of the classical Pascal 
array, whose entries are the well-known binomial coefficients (?) [35]. This feature, which is not put in 
evidence in the traditional description of a Riordan array (see Section [TJ, is highlighted by the umbral 
syntax in the next definition. 


Definition 3.1. Given any umbrae o and 7 , the pair ( 7 , 0 ) denotes an infinite lower triangular matrix 
whose entries ( 7 , a) nz k are given by 


( 13 ) 


(7 , a ) n,k ^ , (7 + k.a) 


i—k 


for n> k > 0 . 


We say that ( 7 , a ) is an exponential Riordan array generated by a and 7 . Exponential Riordan arrays 
of type (p, g), for p, g G f?[A], are defined analogously. Note that if e 7t ~ d(t) and te at ~ h(t), with 
d( 0 ) = 1 , h( 0 ) = 0 , and h'( 0 ) = 1 , then we get a more familiar presentation for (fT3l) : namely 

h(t) k 


r n ,k = (7 , a ) n,k = E 


(7 + k.a) 


i—k 


= n\ [t n } d{t) 


k\ 


In the umbral notation, the identity array is represented by (e,e) and the Pascal array by (v,e), where 
v is the unity umbra and e is the augmentation umbra. Note that, given any umbrae 7 and a, we have 
( 7 , o) njW = 1 for all n > 0. Because of this property, we say that ( 7 , a) is normalized. Any Riordan array 
(d(t), h(t)) with d( 0 ) 7 ^ 0 , h{ 0 ) = 0 , and h'{ 0 ) 7 ^ 0 can be normalized in a natural way; hence restricting 
our attention to normalized Riordan arrays is not a loss of generality. From now on, we will simply say 
Riordan array to mean normalized exponential Riordan array. 

Let /(f) and g(t) be the exponential g.f.’s of the sequences (/*) and (g{), with /o = go = 1- In classical 
terms, a Riordan array ( d(t),h(t )) relates /(f) with g(t) by means of (d(t), h(t))f(t) = d(t)f(h(t)) = g(t). 
This identity is known as the fundamental theorem in the theory of Riordan arrays. In umbral terms, 
we shall write (7,0)77 = w, for 77 and w in d such that f v (t) = f(t ) and / u (f) = g(t). The next result is 
a direct consequence of identity (0) and Definition 13. II 


Theorem 3.1 (Fundamental Theorem of Riordan arrays). For all umbrae 0,7,77 we have 
( 14 ) (7,0)77 = 7 + r}.p.a v . 


The fundamental theorem of Riordan arrays is nothing but a reformulation of the umbral Abel identity. 
Observe also that 

n n 

((7, 0)77)" ~ ^(7, O ) n,k 11 k ^ r ^,kfk ■ 
k =0 k =0 

Hence, after applying the evaluation map E , formula m gives the summation formula 

n 

^2 fk = n! [*"] rf (*)/(*'(0) • 

k =0 
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By taking 77 = v in (fill) . the resulting umbral polynomial 7 + 0.0.^ represents the sequence of numbers 
whose n th term is the sum of all entries occurring in the 77 th row of ( 7 , a); namely, 

n 

(7 + ^Op)" - '^2h,a)n,k ■ 
fc =0 

With the same ease, the product of Riordan arrays is encoded in umbral terms as follows. 

Proposition 3.1. Let (7 ,a) and (a, p) be Riordan arrays. We have 
(15) ( 7 ,a)(cr,p) = (7 + o.0.a^ , a + p.0.a„). 

Proof. Let 77 be any umbra. From the fundamental theorem and identity ©, we have 
( 7 ,o)(cr,p )77 = ( 7 ,a)(o- + 77 ./ 3 .p c ) = (7 + o.0.a v ) + p.0.a + p^-0-a^ 

= (7 + <7.0.a v ) + r).0.{a + p.0.a v ) v = (7 + o.0.a v , a + p^.a^p 
Formula (fl5l) follows from here. □ 

Setting / 7 (f) = d{t), tf a {t) = h(t), f a (t) = a(t ), and tf p (t) = b(t), we recover the traditional formula for 
the product of the Riordan arrays ( d(t),h(t )) and ( a(t),b(t )); that is, 

( d(t),h(t )) ( 1 a(t),b(t )) = (d(t)a(h(t)),b(h(t))) 

Theorem 3.2. Denote by 9tio the set of all Riordan arrays. With respect to matrix multiplication, $Hio 
is a group, where the identity array is (e,e), and the multiplicative inverse of any ( 7 ,ot) £ Olio is given 
by ( 7 ,a ) -1 = (£ 7 ,„,£«)• 

Proof. This is a straightforward verification using identity m- □ 

3.2. The Sheffer group. Closely related to Riordan arrays are Sheffer polynomial sequences. Several 
important classes of orthogonal polynomial systems: Hermite, Charlier, Laguerre, Meixner of the first 
and the second kind, are special instances of Sheffer sequences. These polynomial sequences have been 
deeply studied using finite operator calculus (as shown for instance in Roman’s book W)- A treatment 
of Sheffer sequences under the renewed umbral symbolism can be found in mmm- We briefly review 
in this section the connection between Riordan arrays and Sheffer polynomial sequences. 

A Riordan array ( 7 , a ) defines a polynomial sequence s n (x ) = X)L~o (7) a )n,k x k . By the fundamental 
theorem of Riordan arrays, we can write 

s n (x) ~ (( 7 , 0 ) 0 ;)" ~ (p + x.0.a v ) n , for all n > 0 . 

In a classical fashion, the sequence s„( x) satisfies 1 + X)n>i- S ™( ;C )FT = /■y W namely, it is a Sheffer 

sequence. We can see from the fundamental theorem that each Riordan array determines a Sheffer 
sequence and vice versa. In this sense, we say that (s n (x)) is the Sheffer sequence of ( 7 , 0 :). Sheffer 
sequences are endowed with a binary operation called umbral composition. More precisely, if ( s n (x )) and 
(r„(a;)) are the Sheffer sequences of ( 7 , a) and (tr, p) respectively, then the umbral composition of ( s n (x )) 
and ( t n (x )) is the polynomial sequence ([s r] n (x)) defined by 

n 

[s ( 7 , a )n,k 7k{x) for 71 > 0 . 

fc =0 

Expanding v n (x), we may easily check that the coefficient of x k in [ 5 t] ra (:r) is the (n, fc)-entry of the 
Riordan array (7 ,a)(a,p). Hence, by (fl5l) . the polynomial sequence [sr]„(a:) is the Sheffer sequence of 
(7 + o.0. a,, , a + p.0.a v ). Also, it readily follows that the Sheffer sequence x n of (e,e) is the identity 
element for the umbral composition and that (s„(a;)) is the inverse of (r n (x)) if and only if ( 0 , p) = 
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(£ Ti „, £ q ). Thus, the set ©tjcff of all Sheffer sequences endowed with the umbral composition is a group, 
and moreover, it follows from the previous discussion that ©[jeff and Ttio are isomorphic. 


Subgroup (7, a) 

Appell (7, e) 

Associated (e,a) 
Bell (a, a ) 


(7 ,a)n,k 



Q((t+l).q)’‘- 

Table 2 . 


(7, a) Product 

(-1.7,e) (j,e)(a,e) = (j + a,s) 

(e,£ Q ) (e,a)(e,p) = (e,a + p./3.a v ) 

(£„, £„) (a, a) (a, a) = (a + cr.p.a^^a + <7. / 0.a x> ) 

Some Riordan subgroups 


3.3. Riordan subgroups. Riordan subgroups can be neatly described in umbral terms. A sample of 
them is shown in Table [2] The umbral notation allows for easy verifications. For example, a straightfor¬ 
ward use of m yields (7 ,a) = (7, e)(£, a), showing that any Riordan array is the product of an Appell 
array with an Associated array. Similarly, one can write (7, a) = (7 + (—l).a, e)(a, a), which shows that 
any Riordan array is also the product of an Appell array with a Bell array. Likewise, the fundamental 
theorem of Riordan arrays applied to each of the subgroups shown in Table [2] yields: for the Appell 
subgroup, 

(7, e)r] = 7 + r)./3.e v = 7 + = 7 + r). 

For the Associated subgroup, 

(e, a)r/ = £ + = 77 . 


Finally, for the Bell subgroup, 


(a, a)rj = a + . 


3.4. Recursive properties. Riordan arrays are also characterized by means of recurrence relations. 
The recurrence relations stated next (Theorem 13.31) are direct consequences of Definition 13.11 and the 
umbral Abel identity (HJ). Since these formulas are particular cases of recurrence relations satisfied by 
Riordan arrays in a more general context, we postpone their proofs to Section [5j 


Theorem 3.3. A Riordan array (7,0) satisfies the following recursive properties: 

n—k 


(16) 


(7, a) n , k 


i=0 


n — 1 


(7, a) n -i-i,fc-i , n > k > 1. 


(17) 


(7, a) 




(75 ^)n— l,k—l-\-i 


n > k > 1. 


(18) 


(7 ,Ot) n ,k - 


n—k 

y; (fc.£o«)* (7 j ot) n -k t i , n > k > 1. 
2=0 


Remark 3.1. Traditional formulas are obtained by applying the evaluation map E on the identities m, 
m and (HU). Thus, the umbral equivalence — is substituted by the ordinary equality — and cd, and 
(k.& a Y are replaced with their corresponding moments, say hi ~ a 1 , a* ~ Si l a and a[ k) ~ (&..£<* )*, where 
a( fe) ~ d[f l ](/£ a (f)) fc = *![i*](/a((i/a(i)) <_1> )) fc - The sequence (oj) = (a( x) ) is nothing but the analogue 
of the A-sequence of Rogers [2B] for exponential Riordan arrays. 
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4 . Examples 

The purpose of this section is to analyze concrete examples of classical Riordan arrays and their 
corresponding Slieffer sequences by using the umbral syntax. 


Pascal array. The Pascal array P = (v,e) (written as (e 4 ,i) in classical terms) satisfies 


Pfi , k — 


(■v + k.e) 


n—k 




Since P is an Appell array, we have Pv = v + v = 2, so that, after evaluation, we recover the well-known 
identity 2" = Ylk =o (D• Note that the identity Px = l-\-x encodes the formula (.x + 1)™ = J2k =o (k) xk - 


Stirling arrays. The Stirling array of the second kind, traditionally written as (l,e 4 — 1), is encoded 
by S = (e, — l.t), where t is the Bernoulli umbra (see Table[T]). It is an Associated array. Its entries S n> k 
are the classical Stirling numbers of the second kind, usually denoted by S(n, k). Some known properties 
of this array can be derived from the umbral setting as follows. 

(1) Via generating functions, it is easy to see that (— l.t)® = v. It follows from Theorem 13.11 that 
Sv = £ + v./3.(—l .t) v = /?. After evaluation, the sequence of row sums of this array gives the 
sequence of Bell numbers; that is, o ^( n i = -®n> as expected. 

(2) Let ( 4> n (x )) be the sequence of polynomials in x represented by x./3. Since Sx = x.j3, the Sheffer 

sequence of S is precisely (cj) n (: r)). After evaluation, we obtain <j> n (x) = S( n i &) xk ■ The <f> n (x) 

are called exponential polynomials m- 

(3) Theorem 13.31 allows us to obtain recurrence formulas for the Stirling numbers of second kind. Since 
(—l.t)^ = v, we have (—l.t)^ ~ i(— l.t) 1-1 ~ 1 for all i > 1. Hence, (—l.t)* ~ From identity 
(imil and after evaluation, we obtain (compared with [6]) 

s(n '*) = ?E( n 7 1 ) (rn) s(n - 1 - i '*- 1) - 

Also, since A_i. t = — l.t./3.(— l.i)^' 1> , and (—l.t)^ _1> = u <_1> = y.y, we have A_i. t = —l.t.y. It 
is easy to check that the g.f. of —1 .l.\ is f-i. L . x (t) = Iog ^ +t ) • The rational numbers represented 
by the moments of —l.t.y are known as Cauchy numbers of the first type [Bj. Following m , we 
denote these numbers by ^. From (Ob we get the following recurrence identity that relates Stirling 
numbers of the second kind and Cauchy numbers of the first kind (compare with |19l Theorem 3.1]). 

S(n, k ) = — y] ( ’ . + % S(n - 1 - *, k - 1). 

7-0 V * / 

The moments (fc.J?_i.J* ~ ^ (fe) = d[f*]( log ^_ ) _^ ) fc are known as generalized Cauchy numbers of the 
first type. They show up in the following recurrence relation obtained from (1181) . 


S(n,k ) 


n—k 

J2 c ^ k) S(n-k,i). 

i=0 


Let us now denote by s the array whose entries s Ut k are the Stirling numbers of the first kind s(n, k). In 
traditional terms, we have s = (1, log(l + t)). It is known that this array is the inverse of S. Since S is 
an Associated array and £_i. t = — = t.y, it follows that s = (e, £_i.J = (e, t.x). As shown above, 

known properties involving the numbers s(n, k ) can be obtained through the umbral setting as follows. 
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(1) Via g.f.’s we have [i-x)-D — X-X- By Theorem 13.11 we obtain sv = v.f3.{i.x)i, = v.p.x-X — X- As 
s 0,0 — 1) Si,i — 1 and s Ut o — 0 for n > 1, it then follows after evaluation that 

n 

0 = ^ s{n, k) for n > 2 . 
fc=o 

(2) Note that for any x £ RU A, we have / x . x (t) — (1 +t) E = 1 + J2 n >i( x )™^l’ where ( x) n is the lower 

factorial of x. It follows that ( x.x) n — ( x) n for all n > 0. Since sx = x./3.{— = x.(3.x-X — X -Xi 

the Sheffer sequence of s is the sequence of falling factorials ((x)„), that is 


( x ) n = ^2,s{n,k)x k . 


k =0 

(3) Since (i.x)™ — (x-X)" — (x)n — (—l)" -1 (n — 1)!, identity (fill specializes into 

n—k 


<n,k) = ~Y J 


1=0 


- n (-1) 1 , 

(< + l) 


! s(n — 1 — i, k — 1). 


Also, since M. L . X = i and l 1 ~ bi (the i th Bernoulli number), identity (fTTl) gives 


s(n,k) = -J2 


i—k 


i =0 


k — l + i 


bi s(n — 1 — i, k — 1). 


Finally, for every k > 0, the umbra k.i represents the sequence of generalized Bernoulli numbers 
b[ k) ; that is, {k.i) 1 ~ b (k) = It then follows from (fT51) that 


/ \ n—k 

s(n, k) = b i k) s ( n 


5. Generalized Riordan arrays 

As seen in the Introduction, Riordan arrays can be described in several ways according to the type 
of their generating functions and the manner in which coefficients are extracted from them. Thus, the 
notion of a generalized Riordan array with respect to a given sequence of nonzero numbers arises [35]. By 
means of the urnbral syntax, generalized Riordan arrays are handled in a neat and unified way, extending 
the discussion on exponential Riordan arrays given in Section [3] In this regard, as a direct consequence 
of the urnbral Abel identity (jj), we obtain an important urnbral polynomial equivalence ('Theorem 15.31) 
from which a novel non-recursive formula for Riordan arrays is derived (Theorem [53]). Specializations of 
Theorem 15.41 allows us to recover some known recurrence relations of Riordan arrays (Theorems 15.51 and 
15.61) and unveil a new recurrence relation fTheorem 15.71) . 


Definition 5.1. Let (c„) be a sequence of nonzero numbers such that cq = 1. Let ui be an umbra whose 
moments are uj n ~ )—, for all n > 0. Given two umbrae a and 7, we denote by ^(7, a) the infinite lower 
triangular matrix whose entries ^(7 ,ct) n k satisfy 


(19) 


(7,«) „, fc - 


Cn 

C-k 


(7 + k.a) n ~ k 
{n — k)\ 


for n > k > 0. 


The matrix ^(7,0) is called w-Riordan array and it is said to be generated by a and 7. In general, for 
p,q £ R[A], we define w-Riordan arrays u {p,q) in an analogous way. Observe that w (7, a) n n = 1- This 
means that we are dealing with normalized arrays. 
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Letting e 7 ‘ ~ d(t) = and - W) = J2 n > 1 with d (°) = !> M°) 

we get the traditional presentation 


r n ,k = u {l,ot) nk =E 


Cn (7 + k.a) n k 
c fc ( n—k)\ 


= C n [t n ] d(t) 


h{t) k 

Cfc 


0 , and h'( 0 ) = 1 , 


Exponential Riordan arrays are obtained by setting u> = v (that is, c n = n!), in which case we simply 
write v (j, a) as ( 7 , a). Ordinary Riordan arrays are obtained by setting oj = v (that is, c n = 1), where v 
is the boolean unity umbra whose moments are v 11 ~ n! (see Tabled]). We denote by the set of all 
w-Riordan arrays. All known definitions and properties regarding (exponential) Riordan arrays readily 
extend to w-Riordan arrays. Thus, for instance, it follows that ^Oiio is a group and the natural cu-versions 
of the subgroups of TUo are in turn subgroups of w 9iio. As a sample, we give next a statement whose 
straightforward umbral proof is left to the reader. 


Theorem 5.1. Let u> be an umbra with nonzero moments and let (7 , a) £ Tlio. Then 

-1 -1 

w (7»“) u ( a >P) = J(7,a0(cr,p)] and [ w (7,«)] = j 7 ,a) = £«)• 

Furthermore, the fundamental theorem in the context of w-Riordan arrays is stated as follows. 


Theorem 5.2 (Fundamental Theorem of w-Riordan arrays). For all umbrae 0,7 and 77, we have 
( 20 ) w[Ja, 7)]77 = ( 7 , a)ujr] = 7 + u^./Aoq, . 


If we choose 77 = v in identity (I20[l . we can see that the row sums of ^( 7 , a) are given by 


( 21 ) 


H u (7,<*)„,* - 


k =0 


Cn (7 + 0;./3.Q P ) n 
n\ 


In the same manner as Sheffer sequences are defined from exponential Riordan arrays, we define w-Sheffer 
sequences out of w-Riordan arrays by 


[x = 


/ j OJ 

k =0 


(7 ,a) 


It immediately follows that for every Sheffer sequence (s„(x)), there exists a unique polynomial sequence 
(s%(x)) such that 


( 22 ) oj n s“(x) ~ s n (u>x), n > 0 . 

Since s„(u;x) ~ (7 + ujx.p.a^)™ , it follows from (1^1) that s"(x) ~ n ^ ^We denote by 

71! 

u ©fjeff the set of all polynomial sequences (s%(x)) and refer to it as the w-Sheffer set. 


In order to state recurrence relations for w-Riordan arrays, consider the following simple umbral trick: 
7 + k.a = 7 + (k — m).a + m.a, where k and m are any integers. The following main result is a 
straightforward consequence of the umbral Abel identity ([3]). 


Theorem 5.3. For all n, k £ Z such that n — k > 0, we have 

i—k 


(7 + k.a) n k ~ ^ + (k — m).a + i.\) n k 1 m.a{m.a + (—i).A)* 

z=0 ' 1 ' 


v-' (7 + (k — m).a + i.\) n k 1 to. aim.a + (—*).A)* 1 
- k)\ 2 ^- 

i=0 


(71 — k — i)\ 


i\ 
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Recall that, on the one hand, we have a A ~ m.a(m.a + (—i).A)* 1 , for all i > 1. On the other hand, 
we also have 8. m . aiX = m.a./3. A^ -1> . We then conclude that & m . a , x = It readily follows from 

Definition 15.11 and Theorem 15,31 the next key result. 


Theorem 5.4. 

(23) 


Given any umbra A and any integers m, n, k, with n> k, it holds 

/ \ c n 'C {~! + {k - m).a + i.\) n k 

- L -jj- 

2=0 


(n — k — z)! 


Identity (l23l) is equivalent to 


(24) 


,(7,«)„,* - 


Cfc Cn—k 


( mJi *,*) 1 

2^ Ci -7i- 0,(7 + (* - A )„-fc,i 


2—0 


In particular, the (n, fc)-entry of any exponential Riordan array is given by 


(25) 


(7>«)n,k- L E 


n—k 


i =0 


. — k 


(m.K atX y (7 + (fc — to). a + z.A) 


i—k—i 


Remark 5.1. Theorem 15.41 is important not only for the recurrence relations ('Theorems 15.51 [5.61 and T5.7I) 
that we shall see are derived from it, but also because it provides a way of expressing the (n, £;)-entry of 
any w-Riordan array as a special weighted sum of entries of a related w-Riordan array parametrized by 
A £ A. Formula (l24l) is, to the best of our knowledge, new. Remarkable combinatorial identities as the 
ones shown in Examples 15.11 and 15.21 are direct consequences of this formula. 


Example 5.1. Consider the Stirling array of the second kind S = (e, — l.i) and set A = —l.z in (1^51) . 
Recall that = —l.z.%. Hence, for any n > k > 0 and m G Z, we have 

n—k 


S(n , k) ~ 


E 

2=0 


. — k 


((-m).i.xY ((m -k- i).i) 


i—k—i 


Since ((— m).u.'x) 1 ~ ^ tm) and ((to — k — z).z)” * * ~ (these are respectively the generalized 

Cauchy and Bernoulli numbers that were used in Section [4]), we get the following combinatorial formula 

i—k 


S(n , k) = 


E 

i —0 


n — k 


°n-k-i • 


Now, set A = e in (l24l) . We have = x and X < x> = X> hence = a. The next identity follows. 


Corollary 1. The entries of any w-Riordan array satisfy 

n—k ( \i 

Cn E Ci —jj— J7 + ( k - m )- a ’ £ )n-k,i ■ 


(26) 


;(7,<*)„.* - 


Ck Cn—k n 
2=0 


In traditional terms, formula (|26[) is equivalent to 


7 —k (fc) 

a) 


Tn,k — 


Ck Cn—k . n 
2=0 


E a i 

Ci 7j ^ n—k,i 5 


where (to. a ) 1 ~ a) m) = z![T] (/ a (t)) m and 

ui(Gi a )n,k ~ 

w ( 7 +(fc-TO).a,£:) n _ feii = 


Cn[t n ]f 7 (t) {tfa c Y )) and 

C n -k[t n ~ k ]Mt)f a (t) k - m £. 
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Moreover, setting m = k > 0 in we obtain 

J7,a) n ,fe - 


k (fc.a)* 


C/c C-n—k 


E 

2=0 


( 7 5 £ )n-fe,i- 


•| uj\ I "> ° Jn—k,i 


Example 5.2. Consider again 5 = (e, — l.t) and formula (1261) . For any n > k > 0 and ratZ, we have 




i—k 


2=0 


- — /c 


((— m).i) 1 ((m — fc).t) 


i—k—i I Tl 

.k 


((-kur k ^ry-*i 


In traditional terms, the identity above is written as (see m Proposition 9.1]) 

S ( n,A0=(”)C1. 

By specializing A = a in identity (|24[) and setting m = k > 0, the following linear horizontal recurrence 
relation for w-Riordan arrays is obtained. 

Corollary 2. Any w-Riordan array satisfy the following recurrence relation, 


,(7 ,a). k - 


&k Cn—k 


n—k 

E< 

2=0 




i\ “ 


(7) a ) n ~ 


k.i ’ 


> k > 0. 


Furthermore, by setting A = a in identity (l23l) . Corollary [2] is obtained as a particular case of a more 
general horizontal recurrence relation for w-Riordan arrays. 


Theorem 5.5. For any integer m such that n > m, it holds 


i—k 


J7,«*)„.* - 


Cfc Cji—m 


\' „ (m.R a ) ^ ^ 

Ck-m+i ., ojW) a ) n -m t k-m+i ' 


z—0 


Proof. We can write n — k—i = (n— m ) — (A: — m + *) so that 

c„ (m.Sj)* (j + (k — m + i).a) 

uj\h a )n,k — ~ / . 


(n—m) — (k—m-\-i) 


Ck *7^ *! ((u — to) — (fc — m + z))! 

n—k ( a 

Cn ( to -UqJ Ck-m+i , , 

7T Z—/ 7 wW; a ln-m,fe-m+i ■ 

'-'/c n 6 • Cn—m 

2=0 


□ 


Remark 5.2. Theorem [53] extends identities (fT71) and (1151) of Theorem 13.31 to generalized Riordan arrays. 
Setting m = 1 and applying the the linear functional E in Theorem 15.51 we recover a shifted version 
of the recurrence relation stated in .35} Theorem 5.2]. Moreover, by setting oj = v, we also recover the 
shifted horizontal recurrence relation of m see Theorem 3.1 and Formula (3.3) therein]. More explicitly, 
we have ~ ~ii~’ where a) m) = i\[t l ] (f a ((tf a (t)) < ~ 1) )') m . The sequence Q ,- m) /z! is a generalization 

of the A-sequence of (ordinary) Riordan arrays. 

Now, setting A = e in identity (123|l , we obtain the following vertical recurrence relation for w-Riordan 
arrays. 

Theorem 5.6. For any integer m such that k > m, it holds 

( 7 a ) . ~ Cn Ck ~ m V_ - _ ( m ' a ^ (7 a) 

cd \ 1 1 '-*■ jn,k / v ■ | u) \ I ’ )n—m—i,k—m ' 

^ k ,, ( n—m—2 

2=0 
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Proof. Simply write n — k — i = (n — m ~ i) — (k — m) so that 


;(7,a) 


uUi ^/T i,k 


c n (m.a) 1 (7 + (fc — m).a) 


(n—m—i) — (k—m) 


Ck i! ((n — m — i) — (k — m))! 


Cn Ck—m (^777- • Q?J 1 


Cfc z ~ C-n—m—i 

2=0 


;(T Ot) n —m—i,k—m ' 


□ 


Remark 5.3. Theorem Ei extends identity © of Theorem 13.31 to generalized Riordan arrays. It also 
extends the vertical recurrence relation of ordinary Riordan arrays stated in © see Theorem 3.2 and 
the shifted version (3.3) therein]. Note also that, by setting m = 1 and applying the evaluation map 
E in Theorem 15.61 we obtain a different identity compared with the recurrence relation for generalized 
Riordan arrays stated in [38l Theorem 5.5]. 

As a further example of the power of the umbral symbolic approach, consider A = — l.a in identity 
©. We then obtain the following novel recurrence relation for w-Riordan arrays. 


Theorem 5.7. For any integer m such that 2 k — n > m, it holds 

n—k i a \ i 

(ry a ) Ck—m—i ( f \ 

uj\h )n.k „ / j „ • | uA / ’ /n—m—2i,k—m—i ' 

Ck Cn-m-2i V- 

2=0 

Proof. The equivalence follows by similar considerations as those used in the proofs of Theorem 15.51 and 
Theorem EH noting that we can write n — k — i = (n — m — 2 i) — (k—m —i) and —m.R_ Ua = m.R a _ 1-a . 
The condition 2k — n > m is required to guarantee that k — m — i > 0 and n — m — 2i > 0, in order to 
consider the corresponding lower triangular entries ^(y, a) n _ m _ 2i k _ m _ i , for 0 < i < n — k. □ 


In classical terms, the umbral identity given in Theorem 15.71 is written as 


(27) 


Tl k (m) 

_ Cn \ ^ Ck—m—i ^ 

T'n.k — / . Cn—m—2i.k—m—i •> 

Ck n C n — rn —2i 1- 
2=0 


where ( - m.£_!.„)* ~ a( m) = i\[t l ](f a ( ( t // Q (t)) < 1> )) ? 


Example 5.3. As an ordinary Riordan array, the Pascal array P is written as with c n = 1 for all 

n > 0. Indeed, v v = v and = v.j3.v { ~ 1) = x> and then it follows from (fl2l) that 

((k + l).v) —r ((n + l).x) —-(n + l) n _ fc . 

v 7 n +1 v 7 n + 1 

Thus, we obtain 


(v + k.v) n k ^ ((1 + k).v) n ^ k + 1 (n + l)„_fc _ /n\ 

{n — k)\ ( n—k)\ n +1 (n— k)\ \k) 


Now, the g.f. of the umbra v is fv(t) = so that = t — t 2 and (j^y) < x> = 1 ^ 4t . Therefore, 
we have 

Note that is the generating function of the Catalan numbers. Accordingly, we denote by C- m) 

the generalized i-th Catalan number of order m given by 


Ci m) = 





1 + y/T^t 































A SYMBOLIC TREATMENT OF RIORDAN ARRAYS 


15 


being C^ the classical i-th Catalan number whose explicit formula is (7j = AjO- By setting a^ m) = 

i \ C^ m) and substituting the corresponding values in equation (E71) . we get the formula 

( n \ = V r (m) ( n - m ~ 2i \ 

W U * \ k — m — i ) 

whenever 2 k — n > m. 
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